TILTED LEMAITRE-TOLMAN-BONDI SPACETIMES: HYDRODYNAMIC AND 

THERMODYNAMIC PROPERTIES 



O 

(N 

Or 



O 
I 

W)' 



> 

(N 

o 



X 



L. Herrera *Q A. Di Prisco *| and J. Ibanefl 

Departamento de Fisica Tedrica e Historia de la Ciencia, Universidad del Pais Vasco, Bilbao, Spain 

(Dated: January 25, 2013) 

We consider Lemaitre-Tolman-Bondi spacetimes from the point of view of a tilted observer, 
i.e. one with respect to which the fluid is radially moving. The imperfect fluid and the congruence 
described by its four-velocity, as seen by the tilted observer is studied in detail. It is shown that from 
the point of view of such tilted observer the fluid evolves non-reversibly (i.e. with non-vanishing 
rate of entropy production). The non-geodesic character of the tilted congruence is related to the 
non-vanishing of the divergence of the 4-vector entropy flow. We determine the factor related to the 
existence of energy-density inhomogeneities and describe its evolution, these results are compared 
with those obtained for the non-tilted observer. Finally, we exhibit a peculiar situation where the 
non-tilted congruence might be unstable. 
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I. INTRODUCTION 

It is already a stablished fact that a variety of line el- 
ements may satisfy the Einstein equations for different 
(physically meaningful) stress-energy tensors (see 
[13j and references therein). This ambiguity in the de- 
scription of the source is generally related to the arbri- 
tariness in the choice of the four velocity in terms of which 
the energy-momentum tensor is split. 

Thus, when the two possible interpretations of a given 
spacetime correspond to a boost of one of the observer 
congruence with respect to the other, both the general 
properties of the source and the kinematical properties 
of the congruence would be different. 

This is for example the case of the zero curvature FRW 
model, which represents a perfect fluid solution for ob- 
servers at rest with respect to the timelike congruence 
defined by the eigenvectors of the Ricci tensor, but can 
also be interpreted as the exact solution for a viscous dis- 
sipative fluid as seen by observers moving relative to the 
previously mentioned congruence of observers [4J . An im- 
portant point to mention is that the relative ("tilting") 
velocity between the two congruences may be related to a 
physical phenomenon such as the observed motion of our 
galaxy relative to the microwave background radiation. 

In other words, zero curvature FRW models as de- 
scribed by "tilted" observers will detect an imperfect 
fluid, energy-density inhomogeneity, different evolution 
of the expansion scalar and the shear tensor, among 
other differences, with respect to the "standard" observer 
(see [H for a comprehensive discussion on this example) . 



However, as it has been discussed before (e.g see 
fl6| ) imperfect fluids are not necessarily incompatible 
with reversible processes, accordingly, since the dissipa- 
tive character of the fluid defined by the non-vanishing of 
the divergence of the four-vector entropy flow has an ab- 
solute meaning, it is very important to elucidate whether 
or not entropy production is actually happening. 

In the context of the standard Eckart theory [13] a 
necessary condition for the compatibility of an imperfect 
fluid with vanishing entropy production (in the absence 
of bulk viscosity) is the existence of a conformal Killing 
vector field CKV) \ a such that \ a — wriere V a is the 
four-velocity of the fluid and T denotes the temperature. 
In the context of causal dissipative theories, e.g. [Tc|- 
f23j ], as we shall see below, the existence of such CKV 
is also necessary for an imperfect fluid to be compatible 
with vanishing entropy production. 

Now, as it is well known FRW spacetimes admit CKV 



24]; therefore tilted FRW even though described by im- 
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perfect fluids, might not be associated to irreversible pro- 
cesses (as it seems to be indeed the case, see [HI and [l6[ 
for a detailed discussion on this point). 

It is our goal in this work to study in detail tilted 
Lemaitre-Tolman-Bondi spacetimes (LTB). The reasons 
to undertake such an endeavour are twofold. On the one 
hand LTB dust models [25l - |27j are among the oldest and 
most interesting solutions to Einstein equations. They 
describe spherically symmetric distribution of inhomo- 
geneous non-dissipative dust (see [H, [29[ for a detailed 
description of these spacetimes) . 

They have been used as cosmological models (see [3C 
|34| and references therein), in the study of gravitational 
collapse and the problem of the cosmic censorship [Hj]- 
[42| , and in quantum gravity (43l |44| . 

A renewed interest in LTB has appeared, in relation 
with recent observations of type la supernovae, indicating 
that the expansion of the universe is accelerating. Indeed, 
even if it is true that there is general consensus to invoke 
dark energy as a source of anti-gravity for understand- 
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ing the cosmic acceleration, it is also true that a growing 
number of researchers consider that inhomogeneities can 
account for the observed cosmic acceleration, without in- 
voking dark energy (see [45T452} and references therein). 

On the other hand, we know that LTB does not admit 
CKV [53j and therefore heat flux vector appearing in the 
energy-momentum tensor of the tilted congruence would 
be necessarily associated to a "truly" (i.e. entropy pro- 
ducing) dissipative phenomenon. Indeed, we shall show 
that in the context of causal dissipative theories too, the 
existence of a CKV is necessary for the compatibility of 
imperfect fluids with vanishing entropy production. 

We shall provide a possible explanation of such "truly" 
dissipative processes based on the non-geodesic character 
of the tilted congruence. 

The "inhomogeneity factor" , i.e. the variable repre- 
senting those aspects of the fluid distribution which are 
responsible for the appearance of energy-density inho- 
mogeneities has been identified for the tilted congruence 
and its evolution equation has been integrated. 

Also a discussion about the stability of the non-tilted 
congruence and the consequences of attaining the so 
called "critical point" are presented. 

Finally a summary of the different issues discussed 
throughout the text, is given in the last section. 



II. TILTED LTB SPACETIMES 

We shall now provide a general description of tilted 
LTB spacetimes, as well as the basic equations required 
for our discussion. 

Let us consider a line element of the form 



ds 2 = -dt 2 + B 2 dr 2 + R 2 (d6 2 



% (1) 



where B{r, t) and R{r, t) are functions of their arguments, 
and an energy momentum tensor describing a dust dis- 
tribution with energy density p, in comoving coordinates 



respect to which a fluid element has radial velocity to. 
The corresponding tilted congruence is characterized by 
the vector field 



1 



ui 



(1 -CJ 2 )!/2' B{\-U1 2 ) 1 / 2 



,0,0 



(•5) 



We shall now assume that the source as seen by the 
tilted observer consists of an anisotropic fluid (principal 
stresses unequal) dissipating energy in both, the stream- 
ing out limit (a radially directed flow of a null fluid) 
and the diffusion approximation (described by means of 
a heat flow vector whose four velocity is given by 
(|5|). In this case the energy momentum tensor reads: 

T afi = (fi + P±)V a Vp + P±g a p + (P r - P±)s a sp 

+q a Vp + V a q/3 + el a lp, (6) 

where as usual, fi,P r ,P± denote the energy density, the 
radial pressure and the tangential pressure respectively, 
e is the energy density of the null fluid describing dissipa- 
tion in the streaming out approximation, l a is a null four 
vector and the heat flux vector satisfying q^V^ = 0, is 

(7) 



r 



where 



and 



(1 -W 2 ) 1 ^' B (l_ w 2)l/2 



.0,0 



( 1 + " 1 + " o o] 



(8) 



(9) 



satisfying 

V a V a = -l, V a q a = 0, s a s a = l, s a V a =0, 
l a V a = -l, l a s a = l, l a l a = 0. (10) 

An equivalent form to write the energy momentum tensor 



T, 



(2) 



The general form of LTB metric is obtained from the 
integration of the (t r) component of Einstein equations 
which in turn implies the vanishing of any dissipative flux 
in ([2]), the result of such integration produces 



B(t,r) 



R' 



[1 + k(r)} 



1/2 ■ 



(3) 



where k is an arbitrary function of r and prime denotes 
derivative with respect to r. 

In the above we have assumed the congruence to be 
comoving with the fluid and accordingly 



V>* = (1,0,0,0). 



(4) 



In order to obtain the tilted congruence, let us perform 
a Lorentz boost from the Locally Minkowskian frame as- 
sociated to to the Locally Minkowskian frame with 



T a p = fiV a Vp + Ph a p + U a0 + q (s a V/3 + V a sp) , (11) 
where 

h a p = g a p + V a Vp, (12) 



P = 



P r + 2P ± 



n = p r - p ± , 



n Q/ 3 = n ( s a s/3 - -h af} ) , 



(13) 
(14) 

(15) 



ft, = fi + e; P r = P r + e: q = q + e. (16) 
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A. Relationships between tilted and non— tilted 
variables 



From ([2]) and (fTT|) it is a simple matter to obtain the 
following relationships linking tilted and non-tilted vari- 
ables: 



P r = (J, - (J,, 



1-UJ 2 



(17) 



(18) 



Observe that the expression for the pressure corresponds 
to that of ram pressure, which is intuitively evident. 



2. 0,g = 

In this case it follows from the above equations 
/2 \xuj 



P r = fJU, fl 



1-UJ 



1-uj 2 



(21) 



q = -e- 



flUJ 



= M 



1-U) 



P r — flUJ 
1-UJ ' 



(19) 



It is worth noticing that for the tilted observer dissipa- 
tive fluxes ( in either approximation) should be present 
since q = e = implies u) = 0. 

In the case of bounded configurations the second fun- 
damental form would be discontinuous at the outer 
boundary implying the presence of a thin shell there. 
Some special cases are: 



1. e = 0,q/0 

In this case it follows from the above equations 
2 P- M w 



P r = \IUJ 



1-UJ 2 



1-UJ 2 



(20) 



Pr = 



In this case we have 



[IUJ 

1-uj' 



(MjJ 
1-UJ 2 ' 



B. Einstein equations 



In terms of tilted variables, Einstein's equations 



G a /3 = R a f3 ~ -^Rg a fi — &KT a p, 



take the form 



(22) 



87rT 0l 



8tt 



1-UJ' 



(p, + P r u? + 2<7w) = I 2 



B R\ R 
B + RIR 



D 



' R 



'~B ~R 



(23) 



8ttT 01 



87TB 

~ ~ 2 



{p, + P r )uj + q(l + uj 2 ) 



R' B R' 
~R ~ 



(24) 



8ttT u 



8irB 2 
1-uj 2 



jlui 2 



P r + 2quj 



-B 



' R 



R 
R 



(25) 



where dot denotes derivative with respect to t. 

Since the Einstein tensor is the same for the tilted and 
non-tilted observers we should have 



/i + P r ui 2 + 2qu> 



(fi + P r )u) + q(l +cj 2 ) =0, 
(27) 



- 2 
jJLU) 



P r + 2qu = 0, (28) 

which follow at once from (fT7]) -(fl~9 |) . Also from ([26|l 

Pi. = 0. (29) 

Fluids with vanishing tangential stresses have been con- 
sidered in the past in different contexts 0, EH . 



III. QUANTITIES DEPENDING ON THE 
CONGRUENCE 



Since we are going to compare the physical picture as 
described by two different congruences of observers, it 
should be obvious that quantities depending explicitly 
on the congruence would play a fundamental role in such 
study. We shall consider two different kinds of quantities, 
namely: kinematical quantities and dynamical quanti- 
ties, these latter being defined in terms of Riemann and 
Weyl tensors. 



A. Kinematical quantities 

In the absence of rotations (as is the case here) the 
congruence is described through the three kinematical 
quantities: the four-acceleration, the expansion and the 
shear. The four acceleration a a is given by 



where 



s/T 



a a = V.fiVP = as c 



D 



{1-LU 2 ) B(l-LU 2 ) 



(30) 



(31) 



It is worth noticing that the tilted congruence is no longer 
geodesic, in contrast with the non-tilted one. This fact is 
going to play a relevant role in the physical interpretation 
of results, as we shall see below. 

It is also worth noticing that while a nontrivial solu- 
tion for a — might exist for a specific LTB spacetime, 



producing a specific "velocity" field (uj), the generic sit- 
uation is characterized by a ^= 0. At any rate we were 
unable to find such a solution. 
Next, the expansion 8 = V.^ is 



e 



t 



B R ujuj J ujR' 

B + R + (I - uj 2 ) + B(1-uj 2 ) + ~BR 

(32) " 



Finally, the shear tensor is defined as usually by 



CTaP = V(a;l3) + a ( a Vl3) ~ g@^a^, (33) 

which in this particular case may also be written as 



where 



and 



r 



(34) 



(35) 



f 



VI -uj 2 



B R 



ujR' 



B R (1 -uj 2 ) B{l-u 2 ) BR 

(36) 

As it is evident from the above, for the non-tilted ob- 
server (w = 0) the fluid is geodesic (a = 0) and the 
expansion and shear take the standard form. 



B. Dynamical quantities 

Dynamical congruence dependent quantities are de- 
fined from the Weyl and Riemann tensors. Thus let us 
first introduce the Weyl tensor, which is defined through 
the Riemann tensor P^g„) the Ricci tensor R a p and the 



curvature scalar 1Z, as: 



C 



R 



-R a pS^ - ~R al j.8p 



+\%9afi + h^(S$g afl - g a pSp. (37) 

The electric part of Weyl tensor (in this particular case 
the magnetic part vanishes due to the spherical symme- 
try) is defined by 



E a p = C ali j3 V V^V' J , 



(38) 
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we may also write E a p as: 



(39) 



where 



1 


R 


B 






R 




2 


R 


B 


ti- 


J) 


R 




1 




r r" 




P' 




2B 2 




ff 


+ P 


m 



2*?' (40) 



Next, let us introduce the tensors Y a p, X a p and iT ai g 
which are elements of the orthogonal splitting of the Rie- 
mann tensor and are defined by [5(| [57} 



Y 



a/3 



R aj 0sVV 5 , 



(41) 



X a/3 =* P* 7 ^V 4 = -V al ep R: p p S VV s , (42) 



and 



^ = lv a7 e P R ep p S VW 5 = -l eoiep V 5 R ep p S , (43) 

where rj aiep is the Levi-Civita tensor, e Qcp = T^ 7 77 7aep 

and i?* / 3 7<5 = 2^p7<5^q/3 £P - 

With these definitions and Einstein's equations we find 

Yi = 47T (f + p) fcg - (4*11 - 5) (X - , 

(44) 

X£ = 4tt fcg - (4.11 + £) (s?s a - \h^j , (45) 

Z a p = -±n<te a rf. (46) 

Tensors Y a p and X a( g may be expressed through their 
traces and their trace-free parts, as 



Y a fi = -Y T h a p + Y T f (saSp - -h a fi 



(47) 



Xap = -^X T h a!3 + X T f (s a sp - -h al 3 ) . (4.S ) 



From (|44| and (|45j) it follows at once that 

Y T = 4tt(/2 + 3P), Ftf = 5 - 47rn, (49) 



• In the absence of dissipation, Xtf controls inho- 
mogeneities in the energy density [58| |. 

• Ytf describes the influence of the local anisotropy 
of pressure and density inhomogeneity on the Tol- 
man mass [58| 

• Yt turns out to be proportional to the Tolman 
mass "density" for systems in equilibrium or quasi- 
equilibrium [58| . 

• The evolution of the expansion scalar and the shear 
tensor is fully controlled by Ytf and Yt [H, [H, HI] ■ 

Another interesting congruence dependent quantity is 
the super-Poynting vector P a , which is associated to any 
dissipative flux present in the fluid distribution. From its 
definition 



Pn 



(51) 



and using pt |) -(|ll) ]) and ([I7 ]) -(fl9" )l we have 



P a = 32tt 2 (jl + P r ) q a = -32tt 



2 -2 
CJfl 



(52) 

Evidently for the non-tilted congruence the super- 
Poynting vector vanishes as it should be for LTB de- 
scribed by a non-tilted observer. 



IV. SOME BASIC AUXILIARY EQUATIONS 

For the forthcoming discussion we shall need the ex- 
plicit form of some basic equations, these are: the equa- 
tions of motion , two differential equations relating the 
Weyl tensor with physical variables and the transport 
equation for the heat flow. 



A. Equations of motion 

The two independent components of Bianchi identities 
Tg. a = 0, after some lengthy calculations can be written 

as 



p,* + fie + tf 



~ 2P' /— 
- | ~ <"> ' 777; \ 1 



2w 



0,(53) 



X T = 87T/J, X TF = -£ - 47711. 



(50) 



These scalars which obviously are congruence dependent 
were introduced in (58j and have been shown to play a 
relevant role in the study of self-gravitating systems, in 
particular; 



and 



PJ,+(fi + P r )a+^ [29 + (T-3w(liiJe) t ]+«*=0 ) (54) 
with /t = f aS ° and /* = l a V a . 
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B. Equations for the Weyl tensor 

As mentioned before two differential equations for the 
Weyl tensor will be needed for the discussion below. 
These two equations originally found by Ellis [6(3, [M| are 
here reobtained following the procedure adopted in (62j 
and expressed in terms of A"tf- They are 



(X TF + 4 7 r/i) t = -3X TF (lni?) t 



+12tt<7 



R 



(55) 



and 



X£ F = -3X TF {lnR)* + 47r<7 



47T/iCT 
1 -0J 2 



+4irq 



2Cj 



r'VT^j 2 

BR 



(56) 



where (l55t has been used. 



C. Transport equation 

In the diffusion approximation (e = 0, q — q), we shall 
need a transport equation derived from a causal dissipa- 
tive theory ( e.g. the Miiller-Israel-Stewart second order 
phenomenological theory for dissipative fluids [Hrlllj]). 

Indeed, the Maxwell-Fourier law for radiation flux 
leads to a parabolic equation (diffusion equation) which 
predicts propagation of perturbations with infinite speed 



(see [22j, |23|, [62 



65] and references therein) . This sim- 
ple fact is at the origin of the pathologies 16611 found in the 
approaches of Eckart [llj and Landau |67l ] for relativis- 
ts dissipative processes. To overcome such difficulties, 
various relativistic theories with non- vanishing relaxation 
times have been proposed in the past [l8l - l2ll . m§L l69j . The 
important point is that all these theories provide a heat 
transport equation which is not of Maxwell- Fourier type 
but of Cattaneo type [7(| , leading thereby to a hyperbolic 
equation for the propagation of thermal perturbations. 

A fundamental parameter in these theories is the re- 
laxation time r of the corresponding dissipative process. 



This positive-definite quantity has a distinct physical 
meaning, namely the time taken by the system to return 
spontaneously to the steady state (whether of thermo- 
dynamic equilibrium or not) after it has been suddenly 
removed from it. Therefore, when studying transient 
regimes, i.e., the evolution between two steady-state sit- 
uations, r cannot be neglected. In fact, leaving aside 
that parabolic theories are necessarily non-causal, it is 
obvious that whenever the time scale of the problem un- 
der consideration becomes of the order of (or smaller) 
than the relaxation time, the latter cannot be ignored, 
since neglecting the relaxation time ammounts -in this 
situation- to disregarding the whole problem under con- 
sideration. 

Sometimes in the past it has been argued that dissi- 
pative processes with relaxation times comparable to the 
characteristic time of the system are out of the hydro- 
dynamic regime. However, that argument can be valid 
only if the particles making up the fluid are the same 
ones that transport the heat. But, this is never the case. 
Specifically, for a neutron star, r is of the order of the 
scattering time between electrons (which carry the heat) 
but this fact is not an obstacle (no matter how large the 
mean free path of these electrons may be) to consider 
the neutron star as formed by a Fermi fluid of degener- 
ate neutrons. The same is true for the second sound in 
supcrfluid Helium and solids, and for almost any ordinary 
fluid. In brief, the hydrodynamic regime refers to fluid 
particles that not necessarily (and as a matter of fact, 
almost never) transport the heat. Therefore large relax- 
ation times (large mean free paths of particles involved 
in heat transport) does not imply a departure from the 
hydrodynamic regime (this fact has been streseed before 
[Tl\ , but it is usually overlooked). 

Thus, the transport equation for the heat flux reads 



Th a0 V^qp. 7 +q a = -nh ap (T p + Tap)-^-nT 



a/3 , 



TV? 

kT 2 



P 

where n denotes the thermal conductivity, and T and 
t denote temperature and relaxation time respectively. 
The transport equation has only one independent com- 
ponent which may be written as 



(q+^+qd-^ 2 









K 




+ ujT ^ 









kT 2 



B 
B 

- ' -(—)'' 
T 2 J B \kT 2 ) 



u;(f-w 2 ) B(1-uj 2 ) 
- T 4q(l-^ /2 - 



(58) 



r 



We are now in capacity to analyze some relevant aspects of tilted LTB spacetimes. 
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V. DOES THE TILTED OBSERVER DETECT A 
REAL DISSIPATIVE PROCESS? (AND WHY?) 

We have seen that for the tilted observer the energy 
momentum tensor corresponds to that of an imperfect 
fluid. However we know that non-dissipative (reversible) 
processes within imperfect fluids are not forbidden a pri- 
ori, in the context of the standard irreversible thermody- 
namics (see [I3 - [l6| and references therein). 



What is the situation for a causal thermodynamic the- 
ory such as Israel-Stewart? In order to eilucidate this 
point we shall for simplicity consider in the dissipative 
part of the energy-momentum tensor only the terms as- 
sociated to the heat flux, excluding any shear and bulk 
viscosity term as well as viscous/heat coupling constants. 
Then from the Gibbs equation and Bianchi identities it 
follows (see eq.(42) in [7g|) 



TS a a = -q° 



(59) 



r 



where S a is the entropy four-current and /3i = -Sf. 

Let us first review the situation for the standard irre- 
versible thermodynamics, in this case we have r = and 
(l59l) becomes 



(60) 



which after simple manipulations takes the form 

1. 



QCt 



- rpap p 
y 1 dis~x9a 



(61) 



where C x denotes the Lie derivative with respect to the 
vector field X a = and T%?_ = V a q + V f3 q a . From 
the above is evident that if \ defines a conformal Killing 
vector (CKV), i.e. 



for an arbitrary function ip, then 
SZ = 0. 



(62) 



(63) 



However we know that LTB spacetimes do no admit 
CKV [53J , accordingly at least in the context of the stan- 
dard irreversible thermodynamics, our tilted observer de- 
tects a real dissipative process. 

Let us now consider the situation for the causal ther- 
modynamics. In this latter case we obtain from (15£ 



not _ 1 rpotfi r ' 
D ;a — 2 dis - Lj x9afi 7J 



q 2 V^T 
reT 2 



(64) 



The equation above allows two possibilities for the van- 
ishing of S a a . Either the two terms on the right cancel 
each other or both terms vanish separately. Now, since 
the second term on the right of contains two phe- 
nomenological parameters (re and r) which are absent 
in the first term, it follows that the vanishing of S a a in 
the first case would imply a specific relationship between 
those two parameters. While this situation is possible, it 



would refer to a specific example and is certainly not de- 
scribing a generic scenario. Therefore we shall consider 
the second case, which requires the vanishing of both 
terms simultaneously. However this is not possible in 
our case since LTB as mentioned before does not admit 
CKV, and therefore we conclude that also in the context 
of the Israel-Stewart theory, there is entropy production 
associated with the heat flow vector q a . 

For spacetimes admitting CKV (e.g. FRW) vanishing 
entropy production requires the vanishing of the last term 
in (|64[) , which implies 



1 



(65) 



with g being the metric determinant and = f ^rgprH 



implying in turn the conservation of T ■ At this 

point we do not know what is the physical meaning (if 
any) of such quantity, nor can we understand why its 
conservation is necessary for the reversibility of the pro- 
cess. 

It is worth mentioning that reversible dissipative pro- 
cesses may occur in collisionless plasma, an example of 
which is the well known Landau damping [73j]. In that 
case, the dissipation is related to electrons whose speed 
in the direction of propagation of an electric wave, equals 
the phase speed of the latter. It is not clear to us if there 
is some link between Landau damping and reversible pro- 
cesses satisfying (|52^1 and (|6"5")) 

The remarkable fact that the tilted observer detects 
a real (entropy producing) dissipative process while for 
the non-tilted observer the evolution proceeds adiabati- 
cally, requires a deeper analysis. For doing that we shall 
heavily rely on a discussion presented in [74[ where it is 
shown that forces may be interpreted in terms of colli- 
sional terms appearing in the Boltzmann equations, and 
thereby producing entropy. Basically, what authors of 
74] show is that a specific collisional interaction may be 
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mapped onto an effective force, implying thereby that 
there exists a certain freedom to interpret collisional 
events (producing entropy) in terms of forces (and vicev- 
ersa). 

Now, we have seen that for the tilted observer the 
congruence of is (in general) non-geodesic, leading 
such an observer to conclude that some "force" other 
than gravitation is acting on the fluid. If we interpret 
this "force" as a collisonal term in the Boltzman equa- 
tion we can understand why the tilted observer detects 
a truly dissipative process while for the non-tilted ob- 
server the evolution occurs adiabatically. In this sense 
it could be said that the four-acceleration is producing 
extra entropy. On the other hand we saw that dissi- 
pation is a distinctive characteristic associated to tilted 
observers, accordingly it could also be concluded that 
those dissipative fluxes are at the origin of the observed 
four-acceleration. 



evolution are different. We shall now elaborate on this 
issue in some detail. 

First of all observe that from (|5"5"j) it follows that, as- 
suming the fluid to be regular everywhere, 



# = <=> p) = 



with 



where 



* = X T f - 



ft AR 3 dx' 



R 3 



A = 12nq 



rVT 



R 



(66) 



(67) 



(68) 



VI. THE INHOMOGENEITY FACTOR AND ITS 
EVOLUTION 

As it is well known, the energy-density in the "stan- 
dard" (non-tilted) LTB is inhomogeneous. The same is 
true for the tilted LTB, however in this latter case the 
physical factors related to that inhomogeneity and their 



and a; is a parameter of the curves of the congruence 
defined by s a . We shall refer to ^> as the inhomogeneity 
factor. 

In the non-tilted case (w = q = 0), we have Xtf = 
—£, and therefore the Weyl tensor becomes the inhomo- 
geneity factor, a well known result [75$ . 

Next, we shall use (f51))) to find the evolution of Re- 
placing (1571) into (|56l and integrating we obtain 



J 



J S 47r<7 t + 47rg(w9 



BR 



R 3 ds' - ft AR 3 dx' 



R 3 



(69) 



where s is the parameter of the curves of the congruence 
of V a . In the non-tilted case (|69[) becomes 



£ = - 



47r f Q n<rR 3 dt 
R 3 ' 



(70) 



implying that deviations from an initially homogeneous 
configuration depend on the shear [75] . We see that in the 
tilted version the situation is by far more complicated, 
and deviations from an initially homogeneous configura- 
tion depend also on the dissipative flux. 



Then, immediately after perturbation ("immediately" 
understood in the sense above), 



w = q = 0; 



Cj ~ q 7^ (small). 



(71) 



With the above conditions wc obtain from ([54")) evalu- 
ated just after the perturbation 

fiu + q = 0. (72) 

Next, evaluating the transport equation (|58[) just after 
the perturbation, we obtain 



VII. ON THE STABILITY OF THE 
NON-TILTED CONGRUENCE 

Let us consider a non-tilted congruence, which at t = 
is submitted to perturbations keeping the spherical sym- 
metry. For simplicity we shall consider the possibility of 
dissipation only in the pure diffusion case (e = 0). 

We shall study the perturbed system on a time scale 
which is small as compared to the thermal relaxation 
time and the hydrostatic time scale. 



rq = —kTlj, 



(73) 



where use has been made of the fact that just before 
the perturbation the system is geodesic and in thermal 
equilibrium, and therefore T' = 0. Combining (|72|) with 
([73")) we obtain 

kT 

) = tbfjiCL - a) = 0, (74) 

r 



with a = — . 

fiT 
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From the above it is obvious that if a ^ 1 then lo — 
and taking repeatedly, time derivative of (j54")l and (|S"8|) it 
follows that time derivatives of any order of w vanish, im- 
plying that the non-tilted congruence can be analytically 
extended beyond t = 0. 

However if a = 1, it is not longer possible to assure 
the stability of the non-tilted congruence after perturba- 
tions. The situation described by such a condition has 
been studied in detail in the past, (see [HJ], (76|-[83| and 
references therein for details). Basically that condition 
(usually referred to as the "critical point") implies the 
vanishing of the effective inertial mass density and be- 
cause of the equivalence principle, of the passive gravita- 
tional mass density, leading to important consequences 
in the dynamics of gravitational collapse. 

In order to evaluate the circumstances under which 
such condition appear, observe that in c.g.s. units 



kT=£[k][T\, r 



(75) 



where G is the gravitational constant (G = 6.67 x 
10 -8 <7 _1 cto 3 s~ 2 ) and [k], [T], [t] and [fj] denote the nu- 
merical value of these quantities in erg s" 1 cm~ 



s and g cm 3 respectively. 
Thus 



kT 

TfJ, 



1 MP1 „ 1n -40 



81 



x 10" 



(76) 



At present we may speculate that a may increase 
substantially (for non-negligible values of r) in a pre- 
supernovae event. Indeed, at the last stages of massive 
star evolution, the decreasing of the opacity of the fluid, 
from very high values preventing the propagation of neu- 
trinos (trapping [84^. to smaller values, gives rise to neu- 
trino radiative heat conduction. Under these conditions 
both k and T could be sufficiently large as to imply a 
substantial increase of a. In fact, the values suggested 



in [H ([«] fa 10 37 ;[T] 10 13 ; [t] : 
c.g.s. units and Kelvin) lead to a 



10 
1. 



-4. 



; [iA « io 1 



Next we have isolated the inhomogeneity factor, which 
differs drastically from the corresponding factor in the 
non-tilted case. There too the dissipative fluxes makes 
the difference. This result appears to be relevant with 
respect to the Penrose's proposal [86] to define a grav- 
itational arrow of time. Indeed, since the rationale be- 
hind Penrose's idea is that tidal forces tend to make the 
gravitating fluid more inhomogeneous as the evolution 
proceeds, thereby indicating the sense of time, it should 
be clear that all factors associated to energy-density in- 
homogenity (and not only the Weyl tensor) should be 
present in any definition of the gravitational arrow of 
time a la Penrose, implying thereby that such definition 
would be also congruence dependent. 

Finally we have shown that under extreme conditions 
(the critical point) the non-tilted congruence might be 
unstable, meaning that if that condition is attained, the 
"natural" version of the model would be a non-tilted one. 

We would like to conclude with three remarks: 

• It should be emphasized that our goal here is not 
to provide specific models for given astrophysical 
scenarios, but just to bring out the relevance of the 
role of the observer in the description of physical 
phenomena. 

• Since the physical interpretation of both models 
(tilted and non-tilted) is so different, one could ask 
what interpretation is the better one? However we 
agree with Cooley and Tupper Q, in that the key 
issue is not: what the "correct" interpretation of 
the model is? since both are physically viable. The 
point is that each interpretation is related to a spe- 
cific congruence of observers, and the subjective el- 
ement ensuing from any specific choice brings out 
the relevance of the observer in the description of 
a physical phenomenon. This should not be taken 
as weakness of the theory but quite the opposite as 
expression of its richness. 

• In a recent work (53j some of us tried to generalize 
LTB as to admit dissipative fluxes, here we have 
seen that a simple way to do that is just to look at 
LTB from the point of view of a tilted observer. 



VIII. CONCLUSIONS 



We have seen that LTB spacetimes as seen by a tilted 
observer exhibit physical properties which drastically dif- 
fer from those present in the standard non-tilted LTB. 
Particular attention deserves the occurrence of dissipa- 
tive fluxes which are associated to "real" (irreversible) 
dissipative processes. We put forward a qualitative ex- 
planation for the presence of such processes based on the 
equivalence between forces and collision terms discussed 
in [7H and the fact that the congruence of U M is non- 
geodesic. 
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